Averaging theory for the structure of hydrauhc jumps and separation 

in laminar free-surface flows 
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We present a simple viscous theory of free-surface flows in boundary layers, which can accommodate 
regions of separated flow. In particular this yields the structure of stationary hydraulic jumps, both 
in their circular and linear versions, as well as structures moving with a constant speed. Finally we 
show how the fundamental hydraulic concepts of subcritical and supercritical flow, originating from 
inviscid theory, emerge at intermediate length scales in our model. 

PACS numbers: 47.20.Ky, 47.35. +i, 47.32.Ff, 47.15.Cb 



Despite the classical nature of the subject, the flow 
of a viscous fluid with a free surface presents many un- 
solved theoretical problems, even under laminar condi- 
tions. To a large extent this is due to the lack of approx- 
imate methods for describing flows containing separated 
regions, i.e. regions in which the flow is reversed with re- 
spect to the mean flow. Hydraulic jumps are examples of 
such flows. They are large, sudden deformations in the 
free surface of stationary flows and no theory exists for 
their structure — save the full Navier-Stokes equations 
combined with the free-surface boundary conditions, for 
which even numerical solution poses large problems. In 
this Letter we present a method for determining some 
of these flows, which includes viscosity and variations of 
the velocity proflle. Stationary states are obtained as 
trajectories in a simple two-dimensional phase space. 

An inviscid theory of hydraulic jumps, which is still 
the standard hydraulic approach to the subject, is due 
to Lord Rayleigh in 1914 ||^,^. He regarded hydraulic 
jumps as discontinuities (shocks) which can occur in the 
shallow water equations Across a jump the flow de- 
celerates from a rapid supercritical flow, in which distur- 
bances propagate only down stream, to a subcritical flow, 
in which they propagate in both directions. 

The circular hydraulic jump is easy to study experi- 
mentally and to maintain in a laminar state. Here a jet 
of fluid falls vertically onto a horizontal surface (Fig. 1). 
The fluid spreads in an axisymmetric way, and a hy- 
draulic jump is formed at some rj . The value of rj cannot 
be found by the standard theory, and it depends strongly 
on viscosity i/ [||. Further, experiments show clearly that 
a separation bubble, or a recirculating region, forms on the 
bottom in conjunction with the jump ||^~^. 

Separation per se has been studied more intensely in 
boundary layers close to solid bodies, e.g. airfoils, im- 
mersed in a high Reynolds number flow. This line of 
research descends from Prandtl's seminal work j|] in 
which he introduced the boundary layer approximation 
— a radical simpliflcation of the Navier-Stokes equations. 
Prandtl's equations are valid in a thin layer near solid 



(no-slip) surfaces, where the fluid motion is predomi- 
nantly along the surface. But in general, they become 
singular at separation points |]l0|] , where the assumption 
of forward flow breaks down. The boundary layer equa- 
tions can be further simpliflcd by using an averaging tech- 
nique of von Karman and Pohlhausen in which the 
tangential velocity profile is approximated as a low or- 
der polynomium. This model is useful up to a separation 
point, but solutions beyond the point tend to diverge and 
are discarded. Such troubles can be cured by taking into 
account the feed-back effect from the boundary layer on 
the external flow. The "inverse method" jll] makes it 
possible to calculate flows with separation bubbles, and 
analytical results have been obtained for the structure of 
separation points at large Reynolds numbers ||l2| . 

It is natural to employ the boundary layer approxi- 
mation to describe hydraulic jumps since the fluid moves 
nearly parallel to the bottom surface. To avoid the singu- 
larities near a jump encountered in earlier work [^,^ we 
use the Karman-Pohlhausen method. The standard hy- 
drostatic approximation then provides the link between 
pressure and layer thickness analogous to the feedback 
mechanism in the inverse method used above |p^ . 

For a stationary, radially symmetric flow with a free 
surface the boundary layer equations take the form 



wuz — —gh' 



Ur + u/r + Wz 







(1) 



(2) 



where u{r, z) and w{r, z) are the radial (r) and ver- 
tical (z) velocity components, h[r) is the height, and 
we assume hydrostatic pressure. Surface tension has 
been neglected, since it does not appear to be deci- 
sive in determining the structure of the flow, although 
it is necessary for the stability of the flows as dis- 
cussed below. The boundary conditions are no-slip on 
the bottom: u{r, 0) = w{r, 0) — 0, no stress on the 
top: Uz{r,h) = (strictly valid only for small defor- 
mations \h'\) and the kinematic boundary condition at 
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the top: w{r, h) = u{r, h)h' , which ensures the mass con- 
servation: 2'Kr udz =const.= Q = 27: q. By rescal- 
ing the horizontal and vertical lengths and velocities by 
L = [q^v-^g-^f/^, H = [qvg-^fl\ and V = {qi^g^Y^^ 
respectively, all parameters are eliminated from (|l|,^ . 

We average these equations over z, but in contrast 
to earlier approaches ||,^ we shall not assume a self- 
similar velocity profile. Instead the velocity profile is 
parametrized as 



u{r, z) — v{r) (ai(r)7] -t- a2{r)rf' + az[r)rf 



(3) 



where rj = z/h(r). Implementing the boundary con- 
ditions reduces the parameters to just one: A(r), i.e. 
ai = A -I- 3, 02 = -(5A -f 3)/2 and = 4A/3. Thus 
the profile is parabolic when A = and separation oc- 
curs for A — —3. With these assumptions the averaged 
momentum equation (pi) takes the form p3| 



v{F2{\)vy ^ -hf -{\ + 3)v/h^ 



(4) 



where i^2(A) = {J^ u^dz)/ (hv^) = 6/5 - A/15 + AV105 
and rhv = 1. To determine A, one more equation 
is needed, and this is (as in the standard Karman- 
Pohlhausen approach [|j) taken to be (|l|) evaluated on 
the bottom (z = 0), which gives 



h' = -{5X + 3)v/h^ 



(5) 



When V = 1/ (rh) is inserted into (^,|^) we obtain a 
non-autonomous flow for the two-dimensional vector field 
{h, A). The flow has singularities only on the lines h = 
and X — 7/2. It is thus possible to obtain both separation 
and parabolic profiles without crossing singularities. 

We solve the system with two boundary conditions. 
Since the velocity profiles are not measured, we impose 
two surface points hi{ri) and ft.2(j'2), read from a recent 
measurement of Ellegaard et al. [Q. Iterative adjustment 
of A at one end converges to a solution which passes 
through the two chosen points. Figure ||(a) shows com- 
parison of the calculated height h{r) and the measure- 
ments, for two different /12 values. The surface profiles 
near rj show fair agreements, considering the simplicity 
of the model, but rj is off by around 15%. Figure ||(b) 
shows the calculated A(r). A separation zone (A < —3) 
occurs just behind the jump and its size increases as /i2 
is raised, just as observed. The streamlines and velocity 
profiles are determined from A, and this leads to a graph- 
ical representation of the flow as shown in Fig. |l|(b) . 

The model captures the experimental feature that hir) 
inside the jump is little affected by the change in /i2. 
The curves in Fig. ||(a) apparently follow a single curve 
inside the jump, from which trajectories diverge when r 
is increased. Backward integration from r2 automatically 
settles down to this value of A, which helps us since the 
entrance velocity profile need not be specified. Further 
details of the phase space structure can be found in [O . 



Another quantity measured in is the surface velocity 
U{r) = u(r, h{r)). Comparison is made in Fig. |3[ which 
shows quantitative agreement, although rj again comes 
out smaller. There is no free parameter other than /ii,2, 
taken from the experiments. 

One can apply these methods to time-dependent flows. 
Since the time-dependent circular jumps typically in- 
volve breaking of the radial symmetry |^, we take, as 
an example, the two-dimensional (Cartesian) flow down 
an inclined plane. There exists a large body of liter- 
ature p^-p7[| on such flows to which we shall be able 
to compare. We non-dimensionalize in terms of the 
parabolic laminar solution |]l7|] with a constant height 
ho, mean velocity wq, and flux go related by qo = h^VQ = 
ghg sinQ;/(3j^), where a is the bottom slope. The Reynold 
number is i? = voho/iy — q^/v while the Froude number 
is F = Vo/^gho cos a) — i?tanQ;/3. We obtain ||l^ 



ht + {hv)^ = 



^[[hv), + {hv^F2{X)) 



hx cot a 
(A + 3)w/(3/i') 



/i^cota= l~ {5X + 3)v/i3h^) 



(6) 



(7) 



(8) 



where x is the scaled downward distance along the plane. 

We first study stationary solutions to the equations. 
Then, hv — 1 from (|^), and (|^,||) form an autonomous 
two-dimensional system for (A, /i), that can be easily 
studied on a phase portrait. (This is the Cartesian ver- 
sion of (|,|).) There is a unique fixed point /i — 1 and 
A = 0, and thus one cannot find stationary states con- 
necting two different states with constant (/i, A). (This 
can, however, be done for traveling waves; see below.) On 
the other hand an interesting solution [jisl is represented 
by the stable manifold of the fixed point emerging from 
/i = as shown in Fig. ^. The first part of the trajectory 
has ft-j. nearly constant i.e. h « A{x — Xq)/R and 
A w —3/5 until it suddenly jumps up to the fixed point 
values. Inserting into (^,0), we get A = 2.4/F2(— 3/5) ~ 
1.93. We believe this represents flows that are observed 
behind sluice gates though we treat the flow as laminar 
pO| . The conventional hydraulic theory predicts that 
a jump occurs behind a gate when the bottom slope is 
"mild" (i.e. a being less than a critical slope). Corre- 
spondingly, the jump structure in Fig. ^ disappears as 
R tan a — 3F is increased beyond A . 

It is also possible to find a traveling wave solution |T^ ] 
which connects two parabolic laminar solutions of height 
hi at X — ~oo and /i2 a-t a; = 00 pl[ |. These two limits 
thus carry different fluxes, such that the flux is conserved 
in the moving frame only. By choosing the characteristic 
height appropriately, we may set /ii/i2 = 1 without loss 
of generality. Then we define the moving frame by ^ = 
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X — ct, and look for a stationary solutions in ^, which by 
must satisfy c = /i^ + /ii/i2 + ^2(> 3). There are two 
fixed points h = hi^2, both with A = 0, and a heteroclinic 
solution from hi to /i2(< hi) as ^ increases can be found 
@ iff i? tana < 60/i?/(25c2/if-61c/if + 33). Such river- 
bore like solutions are calculated and shown in Fig. ^ for 
a fixed hi and a and varying R. Note that the velocity 
profile always remains near parabolic as A departs only 
slightly from zero, and that the width of the "shock" is 
much larger than the thickness of the layer, unlike the 
steady jump (Fig. ^). 

Finally, we study the dispersion of small disturbances 
in the time-dependent system. The spectrum of the 
uniform state (ft = w = 1, A = 0) allows the dis- 
tinction between super- and subcritical flows, which is 
fundamental to hydraulics but not obvious for viscous 
flows. Surface tension is necessary for the stability calcu- 
lations. An additional term iA-RWh^xx/'i) thus appears 
on the right hand sides of (0,^, where W = n /{ph^v^) = 
%(7v^ / [pgh^siT?' a) is the Weber number. 

Assuming that all disturbances vary like exjp{ikx—iujt) , 
we obtain [ p^ two dispersion branches. In the /c — > 
limit they behave as 0J+{k) ^ 3fc -f ik^{hR/A — cot a) + 

0{k^) andw_(fc) 14/c/25-12i/5i?-|-0(fc3). Thus, the 

uj- branch moves backwards, and the flow is, irrespective 
of the Froude number, "subcritical". From the imaginary 
parts, both branches are stable for a small R, but the 
branch becomes unstable for Rtana > 4/5. This is in 
qualitative agreement with other models, notably those 
coming from perturbation expansions [ p^ , p^|j2^ and from 
averaging The so-called "Shkadov model" is 

identical to our system (^,£|) with a rigid parabolic profile, 
i.e. A = 0, and omitting (pf). 

For very large k the model shows unphysical behavior, 
as one branch becomes unstable. A priori we have no 
reason to expect our model to be well defined at length 
scales much smaller than the normalized height, since 
our starting point is the boundary layer approximation. 
High-frequency oscillations are expected not to penetrate 
far into the fluid, but our assumption of the hydrostatic 
pressure still connects A and h rigidly. We can remedy 
this by modifying (^), such that A depends on a spatial 
average of h and h^ over an interval of the order of a 
fraction of h [ p3| . In this way, the limit k ^ oo now 
corresponds to the Shkadov model jl^ and is stable. 

For intermediate fc, the dispersion behavior is very in- 
teresting. When i?tana > 20/11 w 1.82, the group ve- 
locity of both branches will have positive real parts, and 
the flow is "supercritical" . In terms of the Froude num- 
ber, this inequality becomes F > 20/33 which is similar 
to the classical criterion of > 1 even though our the- 
ory includes viscosity. For large _Rtana, the subcritical 
range of k is so small that this supercritical behavior 
dominates. In Fig. ^, we show dispersion curves (the real 
parts of w±) for = 25, a = 5 [deg], and W = 0.01. The 
slope of the oj+ branch is always positive. On the other 



hand, the slope of the w_ branch changes its sign. The 
subcritical small-fc region is small, and such long-wave 
disturbances become hard to create. Note that the crite- 
rion i?tana = 3F < A 1.93 for the existence of a sta- 
tionary jump (Fig. ^ is almost equivalent to the demand 
that the final state ft = 1 be subcritical (i?tana < 1.82). 
On the other hand, the linearly increasing part before the 
jump in Fig. ^ is expected to be supercritical, but the dis- 
persion around the solution is hard to obtain due to its 
finite extent and the non-uniform character. In contrast, 
for the moving jumps in Fig. ||, super- or sub-criticality 
must be determined with respect to the jump. It can be 
shown ||l^ that the fiow is supercritical in front of the 
jump and subcritical behind, as expected. 

To conclude, we have presented a simple model of free- 
surface flows which can describe separation and the struc- 
ture of the circular and linear hydraulic jumps. 
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It is presumably possible to replace the loss term in (^^) 
to model turbulent flows, as e.g. in [^. 
It is also possible to study homoclinic orbits [|16 
there surface tension appears crucial. 
The model used in [^ has only one branch of dispersion 
relation and thus no distinction between supercritical and 
subcritical flow. 

The modified equation (^ is 1 — (5A + 3)v/3h^ = 
Jl^^ g{x — s)p{s,t)ds where p{x,t) = hxcota — 

RWhxxx/Z and g{x) = oo exp(— Trgga;^). In the go oo 
limit this reduces to {bf) and kc = 2y/T^gi) signifies the 
cutoff in k, which, in Fig. Hwas fee = 67r. 
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FIG. 1. (a) A circular hydraulic jump is formed when 
a liquid jet falls onto a plate from above (photo: courtesy 
of A. E. Hansen), (b) Surface profile, streamlines, and the 
horizontal velocity profile in a cross section, predicted from 
our model (see text and Fig. |^. Note the difference in the 
scales for the two axes. The profile is nearly parabolic at large 
radius, but is strongly deformed near the jump. The shaded 
area is a separation bubble. 
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FIG. 2. (a) Measured h{r) (dashed curves) vs. our model 
(solid curves). The model uses a shooting method from 
r2 = 30 [mm] toward r\ = 12 [mm] for a fixed inner height 
h\ and two outer heights /12. Fluid: 50% ethylene glycol 
(Q = 27[m£/s], u ^l.&x 10~''[m^/s]). Length and velocity 
scales: L = 28[mm], H = 1.4[mm], and V = 12[cm/s]. (b) 
The model predicts a larger separation zone as h2 increases. 
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FIG. 3. Measured surface velocity U{r) (dots) vs. our 
model (dotted lines). Parameters: 80% ethylene glycol 
(Q = 34[m^/s], u = 14.4 x 10"''[mVs]), L = 25[mm], 
H — 1.7[mm], and V — 16[cm/s]. The jump is located at 
rj ~ 24[mm] (experiment) and 20[mm] (model). The inset is 
an enlargement at large r in a log-log scale. 
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FIG. 4. Stationary solutions to the inclined plane equa- 
tions (^-^) with R — 30 and a — l[deg]. Phase portrait (a) 
has a saddle fixed point at {h,X) — (1,0). Dashed curves are 
nuUclines /i' = or A' = 0. Among trajectories (solid curves), 
one stable manifold (drawn thicker) to the saddle point cor- 
responds to the hydraulic jump solution shown in (b) . 
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FIG. 5. Traveling wave solutions to the inclined plane 
equations, hi — l/h2 — 5/4, a — 2[deg], and R = 3,7,9.4. 
The M-profile stays near parabolic (A = 0). Oscillation starts 
when Titan a is near a critical value. 
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